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Abstract 

Generic signatures of quantum chaos found in realistic shell model calcula- 
tions are compared with thermal statistical equilibrium. We show the similar- 
ity of the informational entropy of individual eigenfunctions in the mean field 
basis to the thermodynamical entropy found from the level density. Mean 
occupation numbers of single-particle orbitals agree with the Fermi-Dirac dis- 
tribution despite the strong nucleon interaction. 
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Chaotic dynamics is one of the most extensively developing subjects in physics. Although 
classical deterministic chaos is well understood a rigorous definition of quantum chaos 
does not exist. The relation to the classical limit is straightforward in problems of the one- 
body type. Quantum billiards is the best known example. Classically, regular or chaotic 
features of motion are determined by the shape of boundaries. The corresponding quantum 
level statistics display local correlations and fiuctuations of Poisson or Wigner type. In 
many-body quantum systems a semiclassical picture arises in the mean field approximation 
when the system is modeled by a gas of independent quasiparticles and symmetry (lack of 
symmetry) of the mean field determines regularity (chaoticity) of single-particle motion. As 
the excitation energy and level density increase, the residual interactions between quasi- 
particles transform the stationary states into exceedingly complicated superpositions of the 
original "simple" configurations. Already at early stages of this process the local level statis- 
tics exhibit 0] features of chaos. The pattern of chaotic signatures mixed with the apparent 
failure of the independent quasiparticle model can be called "many-body chaos" 0. 

In this paper we address the question of the relation between the complicated structure 
of eigenstates and the general principles of statistical mechanics. Having at our disposal 
exact eigenf unctions of a model Fermi system with strong interactions (the nuclear shell 
model we compare their statistical properties with those of the equilibrium thermal 

ensemble. 

The statistical approach implies that the observables are insensitive to the actual micro- 
scopic state of the system. Averaging over the equilibrium ensemble should give the same 
outcome as an expectation value for a typical single stationary wave function at the same 
energy 0. This requires the similarity of the generic wave functions in a given energy region. 
Perfect gases give the simplest example of many-body systems where such properties of sta- 
tionary states are evident. To go beyond the perfect gas, we note that the above description 
also fits the notion of stochastic dynamics. In the classical case the correspondence between 
statistical equilibrium and chaotic trajectories exploring the whole energy surface is taken 
almost for granted by many authors, see for example [0. As for the quantum case, already 



the pioneering paper on compound nucleus by Niels Bohr |TT| contains on equal footing 
elements of both patterns, chaos and thermalization. In discussing properties of chaotic 
wave functions, Percival assumes that all of them " look the same" and cover the entire 



available configuration space. According to Berry [jl3[, in systems with the chaotic classical 
limit as a gas of hard spheres, the eigenfunctions should behave like random superpositions 
of plane waves. This conjecture is in fact equivalent to the micro canonical ensemble and 
leads [jl4| to the standard (Maxwell-Boltzmann, Bose-Einstein or Fermi-Dirac) momentum 



distribution for individual particles. 

One can argue that the gas of hard spheres is a specific case of many-body dynamics 
where the interaction is reduced to exclusion of the inner volume of the spheres. However, 



it was shown long ago by van Hove [|15| that the broader class of gas-like systems displays 
quantum ergodicity: a random initial wave function evolves with time into a state which 
gives the same values of observables as the microcanonical ensemble. The assumption of 
randomness or phase incoherence is similar to Berry's conjecture or even to Boltzmann's 
molecular chaos (Stosszahlansatz) . Self-sustained Fermi systems like nuclei behave, accord- 
ing to Landau-Migdal theory, analogously to the gas of interacting quasiparticles. The 
residual interaction cannot be reduced to rare pairwise collisions and is to be treated on 
exact quantum-mechanical basis. The generalization of the results derived for rigid spheres 
to the strongly interacting case is not known. We address this question by comparing the 
signatures of quantum chaos in the nuclear shell model with the pattern of thermal equilib- 
rium. 

The actual computations were performed for 12 particles in the S(i-shell with the Wilden- 



thal interaction ||7|,p!6| which has been tested earlier by numerous calculations of observables. 



Many-body basis states \k) were constructed with good total angular momentum J, its 
projection M, parity vr and isospin T,T^. Within this basis, the hamiltonian matrix has 
diagonal elements which are dominated by the one-body part and numerically are spread 
over the region from -120 to -60 MeV, and off-diagonal elements with an average value of 
about 0.5 MeV. Eigenvalues Ea for states with J'^T equal to 0"'"0 and 2'^0 (with model space 



dimensions = 839 and = 3273 respectively) clearly showed chaotic level statistics [Q. 
The amplitudes Cj^ of eigenfunctions 

\rT;a) = Y.C^\J''T;k) (1) 



have, for a given Gaussian distribution with zero mean value and variance (C^)^ = 1/N°-. 
The localization length A^" gives a measure of local chaos in the vicinity of energy E^. In 
the extreme chaotic case approaches the space dimension N manifesting total mixing 
and delocalization of eigenfunctions. We found A^" ~ 0.9A^ in the middle of the spectrum 
for realistic interaction and A^" ~ A^ for the degenerate model with no stabilizing influence 



of the mean field. Instead of A^", it is also convenient to use the informational entropy [IT? 



s'' = -Y.(ctfH{ctf] (2) 

as well as moments of the distribution function of the components C^. All such characteris- 
tics show that, as the excitation energy increases, the eigenfunctions become more complex 
and the maximum of complexity is reached in the middle of the spectrum. Our measures of 
complexity are basis dependent. We argued [§ that the mean field basis is preferential for 
such an analysis. 

The same process of stochastization can be described in the thermodynamic language. 
A closed equilibrated system with a sufficiently high number of degrees of freedom is ex- 
cited into an energy interval {E^ E + AE) where the density of states with given values 
of exact integrals of motion (J'^T in our case) is p{E). The average ("thermodynamical") 
characteristics are determined by the statistical weight Q{E) = p{E)AE the exact value 
of the uncertainty AE being not important as far as AE ^ E. One can then define the 
thermodynamic entropy S^^{E) = lnQ{E) and temperature T according to 

dS^_l_ 

dE ~ r ^'^^ 

Of course, this description is basis- independent. 

A system with a finite Hilbert space can be heated until the level density saturates at 
maximum entropy and infinite temperature (3). The local level density p{E) for A^ = 839 



states 0"'"0 is presented as a histogram in Fig. 1 together with a Gaussian fit with the 
centroid at Eq = —90MeV and variance aE = 13 MeV. For such a fit, the temperature (3) 
is T = a^/ [Eq — E) . Similar resuhs are vahd for other J'^T classes and the Gaussian fit 
parameters turn out to be the same. This means that one may speak about thermodynamic 
equilibrium. The right half of the spectrum is associated with decreasing entropy and 
negative temperature. 

The Gaussian rather than semicircle p{E) is expected [|] for a many-body system with 
two-body residual interaction. The transition from semicircle to Gaussian level density 
occurs when many-body forces are introduced lifting the selection rules for interactions 
between configurations. On the other hand, the banded random matrix theory predicts, both 
numerically [|1^ and analytically |2^, the semicircle density for a sufficiently wide band of 
nonzero matrix elements around the main diagonal. The realistic hamiltonian matrix is 
banded in the basis of many-body configurations coupled via two-body forces. But the 
matrix is far from being random since its elements are linear combinations of only few (63 
in the sd-shell) two-body matrix elements. 

To compare the global thermodynamic behavior with the properties of individual eigen- 
functions, we have calculated the evolution of single-particle occupation numbers (the 
isoscalar monopole component of the single-particle density matrix) of the orbitals 
A = (Z, j) along the spectrum of stationary many-body states \a), Eq.(l), 

"-S- = ^ 1 4jmr<^limr | «) • (4) 

niT 

The results are shown on Fig. 2 where the panels a,b and c correspond to 0^0, 2"*"0 and 
9^0 (A^ = 657 states), respectively. Although the states within each J'^T class are or- 
thogonal and apparently have nothing in common, all classes exhibit an identical smooth 
behavior of occupation numbers. This is an additional evidence for an equilibrated system. 

In the center of the spectrum all occupancies = ng/(2j -|- 1) become equal to each 
other the common value being 1/2 for our case of 12 particles in the srf-shell. It suggests 
that one can associate to each eigenstate |a) an effective single-particle "temperature" T^_p 
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defined by the Fermi distribution ffj = {exp[{eij — fi)/T^_p] + T^-p changes smoothly 

with Ea being the same for all complicated wave functions near as it should be for an 
intensive thermodynamic quantity. It becomes infinite simultaneously with the thermody- 
namic temperature (3) when the memory of the initial single-particle energies ex is lost. 

The microscopic mechanism of equilibration can be understood from the analysis of 
fragmentation of projected shell model states \J'^T;k). Applying the recipes of statistical 
spectroscopy of French and Ratcliff one can explain the approximately constant occu- 



pation of the si/2 orbital as well as the smooth evolution of occupation factors for ds/2 and 
(i5/2 orbitals as a function of excitation energy. The conclusion is that the thermodynamics 
of the system is determined mainly by the stabilizing action of the mean field, despite the 
strong mixing of configurations. An artificial reduction by a factor 10 of the diagonal matrix 
elements implies (Fig. 2d) constancy of occupation numbers, i.e. vanishing heat capacity. 
In this case one has an analog of a dense hot gas with a very short quasiparticle lifetime. 

Using the occupancies fj^ of individual orbitals one can calculate the single-particle 
entropy of the quasiparticle gas for each state 

= - T.i^j + 1) [fr, In + (1 - fr,) - m- (5) 

Now we have three, apparently different, entropy-like quantities: thermodynamic entropy 
S^'^(E) ~ \np{E), informational entropy (2) and single-particle entropy S^_p (5), the 
latter two for individual eigenstates. In Fig. 3 we juxtapose the corresponding values of 
exp(S') for different physical situations, I, II and III (columns). Rows a,b and c present the 
energy behavior of S^^, S°' and Sf_p, respectively, for 0"'"0 states. 

The column I of Fig. 3 shows the limit of a stable mean field and relatively weak inter- 
action (the diagonal matrix elements are amplified by a factor 10). The thermodynamical 
entropy la is determined solely by the level density. It displays Gaussian behavior of a pure 
combinatorial nature typical for a slightly imperfect Fermi-gas in a finite number of states. 
It is quite similar to the single-particle picture Ic. The informational entropy lb is very low; 
only in the middle region of enhanced level density does one see some effects of mixing. 
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The opposite case III corresponds to the strong residual interaction (as in Fig. 2d the 
diagonal matrix elements are reduced by a factor 10). Almost all states are totally mixed and 
the informational entropy Illb is near its chaotic maximum of exp( 5*") chaotic = 0.48A^ = 
404 for this class of states. S^_p is also constant (IIIc) on the maximum level corresponding 
to the equiprobable population of orbitals. Within the fluctuations, S"" and S'^_p coincide. 
However, the level density still has the Gaussian shape (Ilia) so that the system has normal 
thermodynamic properties. It means that in the absence of the mean field the response of 
the system to thermal excitation cannot be formulated in terms of quasiparticle degrees of 
freedom. Therefore the informational entropy calculated in the quasiparticle basis becomes 
irrelevant from thermodynamic viewpoint. 

The realistic case of the mean field consistent with the empirical residual interaction 
is shown in column II. When the magnitudes are normalized to each other, all three en- 
tropies become identical within fluctuations. Therefore, using Eq.(3), we can identify the 
temperature scales as well. For example, in Fig. 4 we show the temperature which follows 
from Eq.(3) and Gaussian fit (Fig. 1) for p{E), as compared to T^_p extracted from the 
occupation factors of Fig. 2a-c with the use of the Fermi distribution which allows one to 
determine effective single-particle energies e{lj). From the horizontal line corresponding to 
the si/2 orbital, it follows that the chemical potential /i ~ e(si/2). The rf-orbitals can be 
determined with respect to this level. In the a-scale of Fig. 2, the occupation numbers 
have, except for the edges, a constant slope {dnx/da)^. Therefore one can use the middle 
point of infinite temperature and the Gaussian width aE to find 

ca - = Ap{EQ)a% ^.^^ ^ {dnx/da)o. (6) 

Numerical values 6(^5/2) — p = —3.4 MeV and 6(^3/2) — = 4.7 MeV for the effective single- 
particle energies can be compared to the shell model spin-orbit splitting 7.2 MeV near the 
ground state. 

Thus, we have studied the relation of the individual properties of the eigenstates to 
statistical features of the equilibrium thermal ensemble. For the self-consistent mean field 
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and residual interactions, the thermodynamic entropy defined either via the global level 
density or in terms of single-particle occupation numbers behaves in the same way as the 
informational entropy of generic compound states. It allows one to conclude that (i) onset 
of " many-body" chaos is accompanied by the transition to thermodynamic equilibrium; (ii) 
average equilibrium properties of a heated many-body system with strong interactions can 
still be described in terms of quasiparticles and their effective energies in the appropriate 
mean field (this opens the way for explicit calculation of matrix elements between compound 
states 12^ 



Let us stress the special role of the mean field representation ||2^ both for studying the 
degree of chaoticity of specific wave functions and for the statistical description. With 
the artificially depressed or enhanced diagonal matrix elements, the level density still keeps 
the Gaussian shape so the thermodynamic entropy S**^ is qualitatively the same as in the 
realistic case (Fig. 3a). However, with no mean field (Fig. 3-III) the increase of complexity 
measured by the S*" and the mixing of quasiparticle configurations measured by the Sf_p, 
going together, are different from the total heating measured by the level density and the 
"normal" entropy S'^^. The interaction is too strong and the mixing does not depend on the 
actual level spacing. Almost all wave functions "look the same" regardless of level density. 
It means that, with no stable mean field, the quasiparticle "thermometer" cannot resolve 
the spectral regions with different temperatures. 

Finally we would like to give a more formal argument in favor of the direct correspon- 
dence between chaos and thermalization. The general description of a quantum system with 
non-complete information uses the density matrix V which has, in an arbitrary many-body 



basis \k), matrix elements Vkk' = CkC^, where the amplitudes are averaged over the 
ensemble. If the ensemble is generated by interaction with the environment, possible states 
of the entire system are {k;!/) where u characterizes the states of the environment com- 
patible with the state \k) of the subsystem under study. Then Vkk' = J2uCkuCk'u- The 
corresponding statistical entropy S = —Tr (VlnV) is basis-independent and equals zero for 
pure states of the isolated subsystem. For canonical equilibrium ensembles S coincides with 
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the thermodynamic entropy. Let us consider a gas of quasiparticles in the ensemble gen- 
erated by the residual interaction. This makes sense only after proper separation of global 
smooth dynamics from quasirandom incoherent processes. Such a separation defines the 
optimal basis, namely that of the self-consistent mean field (our "simple" states \k)). 
Complicated stationary states \a) mimic the "total" system (quasiparticles + interaction 
field). The ensemble average of Vkk' = Ck^k' to be taken over several neighboring states 
If the amplitudes are uncorrelated and all neighboring states \a) are similar, only 
diagonal elements of Vkk' survive and we come to the informational entropy (2). 

The authors would like to acknowledge support from the NSF grant 94-03666. V.Z. is 
thankful to P.Cvitanovic and V.Sokolov for discussions. 
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Figure captions 



Figure 1. Level density p{E) for 0+0 states; a histogram is compared with the Gaussian 
fit (dashed hne). 

Figure 2. Single-particle occupation numbers, Eq.(4), vs state number a for states O'^O 
(panel a), 2+0 (panel b), and 9+0 (panel c). Panel d shows occupation numbers for 0+0 
states for diagonal matrix elements reduced by a factor 10. For all panels the three curves 
(sets of points) refer to Si/2, (^3/2 and ^5/2 orbitals, from bottom to top. At the middle of 
the spectrum all curves correspond to half-filled orbitals (occupation numbers 1, 2 and 3 
respectively) . 

Figure 3. Entropy-like quantities plotted as a function of energy for 0+0 states. Columns 
correspond to the diagonal matrix elements multiplied by factors of 10 (I), 1 (II) and 0.1 
(III), the latter case coincides with that of Fig. 2d. Rows a, b and c correspond to ther- 
modynamic entropy, informational entropy, Eq.(2), of individual states, and single-particle 
entropy, Eq.(5), of individual states calculated from the occupation numbers, respectively. 
The values in panels c are in units of 10^. 

Figure 4. Temperature found from the Fermi-Dirac occupation numbers of Fig. 2a, points, 
and calculated from the global fit to the level density of Fig. 1, solid hne. 
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